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We present a detailed calculation of the electron-positron production rate using neutrinos in an 
intense laser field. The computation is done for the process f vee via the neutral channel and 
within the framework of the Standard Model. The production rates are tabulated for a range of 
incoming neutrino energies in an electromagnetic field similar to what can be attained by current 
high-intensity lasers. 
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' I. INTRODUCTION 

o : 

' The introduction of a background electromagnetic field can have very important consequences for neutrino physics. 
While the Standard Model neutrino does not itself couple to an electromagnetic field, the effects of such a field 
on neutrinos are made manifest through their interactions with charged particles, real or virtual. Specifically, many 
investigations have calculated the effects of a constant magnetic field on neutrino propagation and interactions. Recent 
areas of interest for such calculations include: the neutrino self-energy 1, 2], the 7-1^-7 vertex i^nucleon interactions 
^1^,1^, neutrino-photon scattering 0jI1i0jEjOiEi0j and neutrino-electron interactions 1l3l. A very nice review 
of these examples and their astrophysical effects has been prepared by Bhattacharya and Pal fl5!| . 

In this article we examine electron-positron pair production through the process v vee. Though normally 
forbidden, the electromagnetic field alters the final states of the electron-positron pair and frees the interaction to 
take place. Such a process could have very significant astrophysical and cosmological ramifications. Gvozdev et al. 

■ have proposed that this interaction could be important in the analysis of supernovas and neutron star coalescence, 
^sj ' and they have offered its role in magnetars as a possible explanation for gamma ray bursts [l6j |. Such effects are 
T-H : important because of the magnetic field strengths associated with these phenomena. For example, neutron stars and 

' supernovae can have fields in excess of 10^^ G 0|, while a magnetar's field is of order < 10^^ G [iq . 

We, however, turn our focus to examining the production of pairs with neutrinos within more terrestrial fields. 
' We investigate the feasibility of detecting 1/ — s- vee events using ultrashort pulsed lasers. Though not approaching 
CIh the astronomical magnetic fields, today's femtosecond terawatt lasers can produce field strengths on the order of 
ji^i E ~ lO^^V/cm {B ~ 10^ G) [l^l- We begin by presenting a review of the derivation of the Volkov field operator 

■ solution for the electron (positron) in an circularly polarized electromagnetic field (Section (nj. Next, we apply this 
Volkov solution to the calculation of the production rate for the process v — > vee (Section IIII|) . Lastly, we provide a 

IL^ ' tabulation of production rates for various incoming neutrino energies in the laser field and estimate the likelihood of 
detecting such a process (Section Hv)l. 
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II. DERIVATION OF FIELD 

As was pointed out in Section any effect an electromagnetic field has on Standard Model neutrino physics is due 
to interaction with charged particles in the interaction. For our process we incorporate the laser's intense field into 
the problem by using the field operator solution of the electron in the presence of an electromagnetic plane wave. We 
begin by solving the Dirac equation for such a case 

{ih^ - e^ - m^c) il){x) ^ Q . (1) 

We choose the electromagnetic field to be circularly polarized and directed along the 2;-direction such that the 4- vector 
potential A{x) is 




: r:i.: (2) 
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where a is the magnitude of the vector potential and hk is the 4-momentum directed along the ^-direction. The choice 
of circularly polarized light as already been pointed out to be a much simpler approach than assuming linearly or 
ecliptically polarized light jl^, El ■ 

D. M. Volkov solved the Dirac equation for an electron in the presence of an electromagnetic plane wave p3. l23j . 
More recent derivations of the solution are outlined in p3. § 40] and We follow the example of Lifshitz, et al. 
[23 | and find that the Volkov field operator tp{x) is 



■ exp 



lea 



2p-k 



p ■ k 
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(pi sin k ■ X — p2 cos k ■ x) 
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In equation (O we have adopted the notation of |25| where the summation over s is a sum over possible spin states. 
We integrate over the particle's momentum p, and the quantities pi and p2 are the components of that momentum 
along the x and y directions, respectively. We note that the momentum that is to be identified with the Volkov state 
is not p. Instead the Volkov momentum, denoted as q, is defined by 



2p-k 



(4) 



Squaring the state's 4-momentum shows that the mass of the Volkov electron has a dependence on the strength of 
the electromagnetic field 



2 , 2 2/2 

me + e a /c 



(5) 



Note that in the absence of an electromagnetic field a ^ the Volkov field operator, equation (j^J, reduces to the 
free-field solution 



lim ip{x) 

a—*0 



(fp 

(2^ 



^/2E^c " ^j2E^c " J 



(6) 



and the Volkov mass m* reduces to the electron mass me 



III. THE RATE OF PRODUCTION T 



A quantity that we wish calculate for this process is the rate of production F. Physically, the rate of production is 
the probability per unit time for the neutrino to emit an electron-positron pair in the presence of the laser field. We 
begin by finding the probability V for the interaction 



(n/ / ^^^E^^E Kp''^«'^';?'-^e;pe-,.d%.,..)r (7) 



where we have summed over the spin states of the final neutrino v' , electron e, and positron e, and averaged over the 
spin states of the incoming neutrino v. The phase-space integral over the final momentums Pe, Pe, and Pi^' has been 
simplified into the form 

JT [ d% 1 _ f d^p,, 1 f d^pe 1 f (Pp-e 1 

i-i-J {27rr2Ef/c J {2TTf2E,,/cJ {2TTf2EjcJ {2nY2E-Jc- ^' 

The scattering operator S for our process is 

^ = 02 ■ 2fl / rf'^V^e(^)7^f-^ + 2sin20w + M^e(a:)Z^(x) 

nP 22 cos2 6'w sm^ J \ 2 2 J 

X / d^y i^Mi" i\ - \A Mv)Za{y) (9) 
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where tpeix) is the Volkov field operator for the electron (equation ©), ip,y{y) is the free- field operator for the neutrino, 
and Z^{x) is the field operator of the Z boson, a is the fine-structure constant, and 9^f^ is the weak-mixing angle. 

In the free-field case the integration over the space-time variables results in a 4-dimensional (5-function that conserves 
energy and momentum. Though this is true of our integral over the space-time coordinate y in equation (|5J), it is 
not true for the integral over x. The Volkov fields V'e(a^) in the J d'^x integral have a more complicated space-time 
dependence. 

We begin by expanding the field operators and integrating out the y dependence. The scattering matrix can be 
written as 



s = -— 
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22 cos2 6lw sin^ 6lw 
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X I (TX \ UsAPe) 1 + 
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X exp 



exp 



-Pel 
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ea 



2pe ■ k 
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ea 



ea 
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(10) 



To fashion equation H1U|I into an analytically integrable form, we employ Hansen's definition of the Besscl function 



g(./2)(t-i/t) ^ ^ j„(z)r. (11) 

n— — oo 

We can recast the left-hand side of equation (|llf) in a fashion more apphcable to the exponentials found in equation (|lUfl 



-iQ sin(fc-a; — (5) _ 
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Since the vector potential A in the 5-matrix has components that go as sin(A: ■ x) and cos(fc • x), we will also need 
to perform a similar expansion for terms like 

sin(fc • x) "j 
cos(fc • x) 

sm(k - x) cos(fc • x) 
cos2(fc • x) 

All of these terms, when expanded, result in forms similar to equation l|12|l . The space-time dependence is factored 
out into an exponential that goes as e"*"*^'^. As an example. 



sin(fc • x)e' 



s'ln^k-x — S) 



C30 

(e^^"+'^'Jn+i(C) - e^("-^^'j„-i(C)) 



(14) 



After these expansions are made, the space-time integral in the ^-matrix reduces to the normal 4-dimensional 
(5-function 



S oc d x exp 



- (qe + qe+ Pv' - Pv - nk) - X 



S oc 5'^{qe + qe+ Piy' - Pu - nk) . 



(15) 
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Physically, equation H15() suggests that it is the Volkov states that are to be considered when conserving momentum 
and energy. Equation IjlSfl also suggests that we interpret the nk term as a number of photons that are pulled from 
the laser field in order to drive this reaction. These photons are what allow this process to overcome the energy and 
momentum constraints that originally prohibited this reaction in the free-field case. In fact, the (5-function can be 
used to apply a constraint on the possible values for n 



2np^ ■ k > 4{mlcf . (16) 

In the proceeding calculation we will take the neutrino and laser beams to be oppositely directed. Therefore, in the 
massless limit of the neutrino the constraint equation (|16|l becomes 



where E^ is the energy carried per photon. We note that the reaction cannot emit photons (n 0), nor can it proceed 
without absorbing some number of photons {n ^ 0). 
The S'-matrix can now be rewritten as 

oo 

S = t{27T)^ S^{qe + qe+P.' -p.-nk) Mn (18) 

n— — cxj 

where we have absorbed the couplings, spinors, Bcsscl functions, etc., into the scattering matrix Mn- For a definition 
of Mn refer to the Appendix. 

Of importance, however, is the square of the S'-matrix. We begin by substituting the S'-matrix (equation (|18|) ') into 
the probability (equation (TJ) 



d^Pf 1 \ f (Ppi, 1 



{2nf2Ef/cl J {2iTf2E^/c 
1 °° 

0^ Yl Yl {'^T^f5^{<le+qe+Pu'-p,.-nk)5^{qe+qe+P,y-Pu~nik)M*^Mn- (19) 



X 

2 

Si/ Se .Se ,s„/ n,m— — CXD 



By inspection, one can see that the two 4-dimensional (5-functions imply that for there to be any contribution to the 
summation, either there is no incoming photon energy [E^ — 0) or, more appropriately, that m — n. Therefore, we 

I 1 2 

can replace A^*jA^„ with the square of the norm of the scattering amplitude \ A4n\ and eliminate the sum over m 



To simplify the result, we have defined the square of the scattering matrix after summing and averaging over spins to 
be 



1-^-1' = E l-^-l'- (21) 
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Next, we can use three of the (5-functions to eliminate the integral over the momentum of the incoming neutrino 



-p- n 



(27r)3 2Ef/c I 2E^/c 



^ {2T:f5\qe+qe+Pi.'-Pu-nk)5{{Ee+Es + E^,~E^~nE)/c) \Mn\ 



X 

n— — oo 
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1 
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This still leaves us with two identical energy-conserving (5-functions. To simplify we follow the solution pointed out 
by and outlined in 3| • We write one of these (5- functions as an integral over time 



oc ^ (^S{{E, + Es + E,, -E,-nE)/c) 

n— — oo 

V y SUE, + Es + E,, ~E,- nE)/c) lim / dt ^^^E,^E,^E^,-E^-r.E)t 

n= — oo I 

and use the other (5-function to reduce the integrand to unity 

oo ^ T/2 

P oc V 8{[E, + E,. + E^, -E^~ nE)/c) — hm / dt 

n=-oo J -1/2 

oo 

P oc y SUEe + E, + E,y -E^- nE)/c) \im T . (23) 

n— — oo 

By making the substitution of equation H23|) into the probability in equation H22|l 

\ f I n— — oc 



we can solve for the quantity of interest in this problem, the rate of production T 

r = 

dt ^ T 

It may be most appropriate to think of the total production rate T as the sum of the production rates for all of the 
individual processes involving n photons 

oo 

r= J2 r„, (25) 



where 



04^) 4^ ^^('^e + ,.+P.-..-n.) \M.\\ (26) 



IV. RESULTS 



Our calculation proceeds by computing the individual production rates r„ (equation (|26|l ) and then summing to find 
the total rate F (equation (|25|l 'l. Though the constraint in equation p7|) imposes a lower bound on the summation over 
n for the total rate, there is no such upper bound. This means that the calculation should include computations of the 
individual production rates for n out to infinity. Fortunately, the individual production rates F„ fall off exponentially 




for large n. How fast (or slowly) these rates fall off depends on what choices are made for the incoming photon 
energy E-y and the magnitude of the vector potential a. For the cases in which the individual rates fall off slowly, we 
calculate the rates to sufficiently large n such that we can characterize the nature of the exponential decay. Using 
this characterization to fit the rates at large n allows us to approximate the summation. 

We choose initial conditions that closely approximate today's high-intensity laser system. These conditions are 
based on the Terawatt High-intensity Optical Research (THOR) laser at the High Intensity Laser Science Group 
at the University of Texas at Austin. The THOR laser is a 20 TW laser beam centered around a wavelength of 
800 nm with an electric field strength \E\ on the order of 10^*^ ~ 10^^ V/cm For our calculations we choose an 
electric field strength of = 5 • 10^*^ V/cm. The choice of wavelength and electric field strength uniquely specify the 
electromagnetic field parameters in our problem. Using the customary relationships, the photon energy is 

Ej = he/ A , 

and the magnitude of the vector potential is 

a = \E\h/E^ . 

The profile of the individual production rates r„ as a function of the number of photons n is given in Fig. ^ for 
incoming neutrino energies of 3 GeV, 30 GeV, and 3 TeV. Notice that for a given choice of photon energy and vector 
potential magnitude, the profiles shown in Fig. Q all follow the same exponential decay. 

The individual rates that make up these profiles are summed in order to find the total production rate (equa- 
tion H25|l 'l. We present the results of this summation in Table for a range of neutrino energies from 1 GeV to 
300 TeV. The total rate of production is also presented for a laser field that is able to generate the same field strengths 
but with a much smaller wavelength of light A = 100 nm. We note that the difference in wavelengths has the most 
significant effect at low incoming neutrino energies. 

Rather than considering the total rate of production for this process, it may have more physical significance if we 
consider the production length A. The production length is the distance that a neutrino must travel in the laser 
field such that its likelihood of producing an electron-positron pair is f — 1/e ~ 63%. That is, the probability for 
production is 

where £ is the distance traveled in the field, and the production length is simply A = c/F. The production lengths for 
the rates tabulated in Tableware shown in Fig. |21 

From inspection of Fig. [3 one can see how feeble this interaction is. For the production length to drop below the 
Hubble length, an estimate for the size of the universe, the incoming neutrino energy must exceed 10 GeV. And at 
energies near 1 PeV the production length is still on the order of a light year. 
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TABLE I: The total production rate V for a given incoming neutrino energy Ei, in the presence of an laser field of electric field 
strength \E\ — 5 ■ lOi" V/cm and wavelength A. 

Production Length 




1 GcV 10 GeV 100 GeV 1 TeV 10 TeV 100 TeV 1 PeV 

FIG. 2: The production lengths A as a function of incoming neutrino energy for a laser field of wavelengths 100 nm and 
800 nm at an electric field strength of = 5 • lOi" V/cm. The dark lines correspond to the Hubble length c/Hq ~ 1.2 • lO'^^ m 
and one light year 1 ly = 9.46 • 10^^ m. 

To further stress the small likelihood of pair production under these conditions, we can consider this interaction at 
a neutrino source such as the Neutrinos at the Main Injector (NuMI) facility at Fermilab. If we assume the most ideal 
set of specifications for the neutrino beam and simplify the beam to be uniform, we estimate the neutrino flux to 
be on the order of 10^^/m^/s per pulse with a pulse every 2 s. The energy profile of the beam is nontrivial. There is a 
low-energy component peaked at 3 GeV, a medium-energy component peaked at 7 GeV, and a high-energy component 
peaked at 15 GeV, but the beam will contain neutrinos with energy up to 120 GeV l^^. However, for estimation 
purposes, we simply assume that all of the neutrinos are at 15 GeV and are uniformly distributed throughout the 
pulse. Using these conditions for the NuMI beam, realistic conditions for the laser beam, and assuming a detection 
region about 1 m long, we estimate that there would be an event once every 10^^ yr. 

Given that our estimate of the time scale for detection of events is twelve orders of magnitude greater than the 
age of the universe, future work must turn towards astronomical and cosmological sources of neutrinos and fields. 
In astronomical phenomena where the magnetic field strengths can reach very high magnitudes, there exists the 
possibility of significant stimulation of this interaction. 
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APPENDIX: THE SCATTERING AMPLITUDE Mr, 

As pointed out in the text, the scattering amplitudes A^„ are related to the S'-matrix by equation 1)18(1 

oc 

S = t{2TT)'^ 6^{qe+qe+Pu' -Pu-nk) Mn- 

n— — oo 

The scattering amplitudes are determined to be 



22 cos2 sin^ 0w V2 2'/ {p, - p,')^ ~ (mzc)^ + iT zmzc^ 

(A.l) 

The C„(C,<5)'s are coefficients that depend on the Bessel functions 

ConiCS) - J„(C)e"' (A.2a) 

cues) = i(j„+i(C)e^^"+^)' + J„-i(C)e^^"-^)') (A.2b) 

CyJCS) = ^(j„+i(C)e'("+i)^- J„-i(C)e'("-i)*) (A.2c) 

C.UC, S) = + \ (jn+2(C)e'("+'^' + J„-2(C)e'("-')*) (A.2d) 

C^y^iCS) = ^ (j„+2(C)e'("+'>' - ^-2(C)e'("-')*') (A.2e) 

C,,„(C, ^5) = i J„(C)e"^ - \ (j„+2(C)e'^"+'^' + J«-2(C)e^^"-''') • (A.2f) 
Recall that C and 5 are defined through equations (|13|l . The matrices are given by 

7^ ( -i + 2sin2% + i7M (A.3a) 



2 2 



2pe-k " "2pe-k 
2pe • fc ° 2pg • k 



= ^^r(;^ (A.3d) 



" 2pe-k °2p,-k^ 2p, ■ k " 2p, ■ k ^ ' 



_ ewfl eafl 



where 7^ and 7^ are Dirac matrices. 
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